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Using a data sample of 6.8 pb' 1 collected with the CMD-3 detector 
at the VEPP-2000 e + e _ collider we select about 2700 events of the process 
e + e _ — y pp and measure its cross section at 12 energy points with about 6% 
systematic uncertainty. From the angular distribution of produced nucleons 
we obtain the ratio Ge/Gm- 


1. INTRODUCTION 


The Born cross section of the process e + e 


by 


—> pp shown in Fig. |T]is given 
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where i/s = 2F , b ea m = -F c .m. is the center-of-mass energy, M p is the proton 
mass, and /3 = yj 1 — 4M//s. The Sommerfeld-Gamov-Sakharov factor [1] 
C = y/(l — e~ y ), y = 7ra/(3, takes into account the Coulomb final state 
interaction. The cross section depends on the electric ( Ge ) and magnetic 
( Gm ) form factors, which are equal at the threshold. To compare different 
experiments, the effective form factor 
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is usually defined. 

I 11 early experiments at the electron-positron colliders [2|, [3|, |4j. |5| in the 
energy range between the proton-antiproton threshold and E cm = 2 GeV, 
the cross section of the process e + e~ —> pp has been measured at six energy 
points only. The accuracy of these measurements is about 25-30% and the 
\Ge\ = \Gm\ assumption is made. In the PS170 experiment at LEAR j(3| the 
measurement of the effective proton form factor and the first measurement of 
the |Ge/Gm| ratio have been performed in the process pp —> e + e~ with 30% 
accuracy. The most accurate measurements of the e + e~ —>■ pp cross section, 
the effective form factor, and the \Ge/Gm\ ratio have been performed with 
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Figure 1: The Feynman diagram of the process e + e —> pp . 


the BaBar [7] detector using the initial-state radiation (ISR) method. How¬ 
ever, PS170 and BaBar results contradict to each other, and new experiments 
are obviously required. 

Additional interest to this energy range is related to an unusual behavior 
of the e + e _ —> 3(7r + 7r _ ) cross section [§J near the proton-antiproton thresh¬ 
old. 

2. THE CMD-3 DETECTOR 


The CMD-3 detector 


0, 0 i 


is installed in one of the two interaction 


regions at the electron-positron collider VEPP-2000 [ill]. The design lumi¬ 
nosity of the VEPP-2000 is 10 32 cm -2 s -1 at the maximum center-of-mass 
energy E c . m . = 2 GeV. The detector tracking system consists of the cylindri¬ 
cal drift chamber (DC) and donble-layer cylindrical multiwire proportional 
Z-chamber, both used for a trigger, and both installed inside a thin (0.2 X 0 ) 
superconducting solenoid with 1.3 T held. The beam pipe inside the DC is 
made of 0.5 mm aluminum with 17 mm inner radius. An inner shell of DC 
is made of a carbon-hber-rcinforced polymer (CFRP) and has 20 mm radius. 
The DC contains 1218 hexagonal cells and allows to measure charged par¬ 
ticle momentum with 1.5-4.5% accuracy in the 100-1000 MeV/c range, and 
provides the measurement of the polar ( 6 ) and azimuthal (</>) angles with 
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Figure 2: Example of an e + e~ — > pp 
event at Eb ea m = 945 MeV. The an¬ 
tiproton stops and annihilates in the 
beam pipe with production of several 
secondary particles. 



Figure 3: Example of an e + e" —> pp 
event at Eb ea m = 970 MeV. The pro¬ 
ton is absorbed in the Z-cliamber and 
the antiproton annihilates with produc¬ 
tion of several secondary particles three 
of which come back to the DC. 


20 mrad and 3.5-8.0 mrad accuracy, respectively. An amplitude information 
from the DC wires is used to measure ionization losses of charged particles 
with (JdE/dx =11-14% accuracy. Two electromagnetic calorimeters (a liquid 
xenon (LXe) one with 5.4 Xq and Csl crystals with 8.1 Ao) are placed in 
the barrel outside the solenoid. BGO crystals with 13.4 A" 0 are used as 
the end-cap calorimeters. The return yoke of the detector is surrounded by 
scintillation counters, which are used to veto cosmic events. 

We use the data samples of 2011 (1.0 T held) and 2012 (1.3 T held) 
runs, collected at twelve beam energy points for an integrated luminosity 
of 6.8 pb~h To study the detector response to a proton-antiproton pair 
and determine the detection efficiency, we have developed a Monte Carlo 
(MC) simulation of our detector based on the GEANT4 [12] package, and 
all simulated events pass the reconstruction and selection procedures. The 
MC simulation includes soft photon radiation by initial electron or positron, 
calculated according to Ref. 13]. 
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3. EVENT SELECTION 


In the studied energy range nucleons have low velocity and high 
ionization losses. According to MC simulation, when the beam energy is 
less than 950 MeV, all protons and antiprotons stop in the beam pipe or in 
the inner DC shell, and antiprotons annihilate producing several secondary 
particles. An example of such an event is shown in Fig. [21 When A c . m . 
is above 952 MeV, almost all nucleons reach the DC sensitive volume, and 
stop in the DC outer shell or in the Z-chamber. Figure [3] shows an example 
of an event, when both proton and antiproton are detected in the DC. For 
beam energies between 950 and 952 MeV, only part of nucleons penetrate 
into the DC volume. According to these differences in the nucleon path to 
annihilation, we use two different approaches for the signal selection. 

3.1. Nucleons reach the DC sensitive volume 

The selection criteria for this class of events are the following: 

a) There are two opposite-charge tracks with the number of DC hits 
A r hit > 4. They are collinear {(56 < 0.25 rad & 5</> < 0.15 rad) or (59 < 0.4 
rad & 5(j) < 0.5 rad) for Ab eam <955 MeV); For Ab eam <955 MeV we have 
very soft momentum of PP particles with high multiple scattering. 

b) The tracks are originating from the beam interaction region within 10 
cm along the beam axis and within 1 cm in the transverse direction. 

c) Momenta of both tracks are close to each other \p± — p 2 \/\pi +P 2 I < 
0.15 (< 0.5 for Abeam <955 MeV); 

d) The total energy deposition in the calorimeters is more than 200 MeV. 

Figure [4] shows a scatter plot of the ionization losses ( dE/dx ) in DC vs 

momentum for a selected pair of tracks. A signal from pp events is clearly 
seen. We require both tracks to have ionization losses above a value, which 
is calculated by taking into account the average dE/dx value and dE/dx 
resolution at the measured momentum. The line in Fig. [4] shows the applied 
selection. 

The distribution of the average absolute value of the nucleon momentum 
for selected pp events at Ab eam =970 MeV is shown in Fig. [5j The number of 
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Figure 4: The dE/dx vs momentum distribution for tracks at Ebeam=970 MeV. The line 
shows the applied selection. 



Figure 5: Average of absolute momen¬ 
tum values of proton and antiproton for 
data at Eb ea m=970 MeV (histogram). 
The line shows a fit described in the text. 



Figure 6: The detection efficiency of 
collinear pp pairs vs £? c . m . for the 2011 
(circles) and 2012 run (squares). 
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background events in the signal range is expected to be negligible. The left 
tail of the distribution is mainly due to the initial state radiation resulting 
in nucleons with a smaller momentum and higher dE/dx value. 

The “declared” collider beam energy was not very precise, and for a few 
energy points the energy was continuously monitored during data taking, 
using the Back-Scattering-Laser-Light system 14]. Based on these measure¬ 
ments and comparing the average momentum value for data and simulation, 
we can determine the c.m. energy with better accuracy. 

We simulate the process e + e” —y pp at the “declared” collider beam ener¬ 
gies, apply the above selections, and fit the average momentum distribution 
of the proton-antiproton pairs with a sum of two Gaussian functions. The ex¬ 
perimental distribution at each energy point is fitted with the corresponding 
MC-simulated function, convolved with an additional normal distribution, 
which takes into account a data-MC difference in the detector resolution. 
The number of events, variance of the additional normal distribution, and 
the momentum difference between simulation and experiment are floating. 
The fit curve is demonstrated in Fig. ED by the line. Using the obtained mo¬ 
mentum difference and proton mass value we calculate a beam energy shift 
-^beam) listed in Table |T] for each energy point. 

The detection efficiency e co n for this class of events is calculated as the 
ratio of the number of selected pp pairs to that of all MC-simulated events. 
The energy dependence of e co u is shown in Fig. ED for two experimental runs. 

To estimate a data-MC difference in the detection efficiency, we select 
a pure class of events with a detected antiproton and check how often we 
reconstruct the opposite proton. We use the following selection criteria: 

- one or two tracks coming from the beam interaction region within 10 
cm along the beam axis and within 1 cm in the transverse direction. 

- one of these tracks has negative charge, has the number of hits Nut >9 
with high ionization losses in DC > 4L™ can _ cr ^), an( j associated energy 
deposition in the calorimeters is more than 100 MeV; 

- the total energy deposition in the calorimeters is from 300 to 1100 MeV. 

Using these selections we obtain the proton detection efficiency for data 
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Figure 7: The correction due to the data- Figure 8: The radiative corrections vs c.m. 
MC difference in the proton/antiproton de- energy, 
tection efficiency. 


and MC-simulation e™, which are calculated as a ratio of the number 
of events with found protons to that of all events with antiprotons. Because 
of the large proton background from beam-gas interactions we cannot select 
a pure sample of detected protons, and assume equal detection efficiencies 
for protons and antiprotons. The squared ratio of efficiencies found for data 
and MC-simulation is shown in Fig. [7] vs c.m. energy, and gives an estimate 
of data-MC difference. We fit these points with a constant and obtain the 
value R = 1.030 ± 0.014 close to unity and it is used as an estimate of the 
systematic error on the detection efficiency. At the point E c m _ = 1901.6 MeV 
we use the correction as determined above, because the effect of material and 
thickness uncertainties of the beam pipe is too big. 

At each energy point the e + e _ —y pp cross section for this class of events 
is calculated from 

N - 

_ iV pP /Q\ 

O'Born — 77 777? > 

Lj ^coll(1 Ojii 


where L is the integrated luminosity, and e co n is the detection efficiency. 
Figure El shows energy dependence of the radiative correction (1-5) calculated 
according to Ref. 13]. 

The c.m. energy , beam energy shift, luminosity, number of selected 
e + e _ —> pp events, detection efficiency, radiative correction, and cross section 
are listed in Table CD 

























3.2. Nucleons are absorbed in the beam pipe or the DC inner shell 

When an antiproton stops in the material of the beam pipe or in the 
DC inner shell, it annihilates with production of several secondary particles, 
which are mostly pions. Part of the produced negative pions are captured 
by nuclei and induce nucleus fragmentation and production of protons and 
neutrons as well as deuterons and tritons. 

Candidates to this class of events are selected with the following criteria: 

a) an event has a vertex with 4 or more tracks located in front of the 
beam pipe or the DC inner shell; 

b) an event has no tracks with energy deposition in calorimeters higher 
than 400 MeV. 

We verify these criteria in a special run without beams, and conclude that 
cosmic events are completely rejected by these criteria. 

The main remaining background is due to the interactions of the particles 
lost from the beams with the detector material. Several pions, protons and 
heavier particles are produced in such interactions. We study this background 
in a special run with one electron or positron beam only, and using data from 
c.m. energy points below the threshold. 

To obtain the number of e + e~ —» pp events we use additional information 
from the calorimeters. The distribution of the energy deposition from the 
background events, shown in Fig. [Deleft) with a fit function, is obtained 
from the runs where the c.m. energy was below the pp production threshold. 
Figure^right) shows the combined distribution of the total energy deposition 
for three c.m. energies above threshold; 945, 950 MeV(run 2011) and 950 
MeV(run 2012). A signal from the antiproton annihilation is clearly seen. 
We describe this signal with an additional Gaussian function and use its 
parameters to obtain the number of e + e _ —> pp events at each c.m. energy 
point. The obtained numbers of signal events vs c.m. energy are shown in 
Fig. HO) 

To obtain the cross section we need to determine the detection efficiency 
for this class of events. 

To calculate the efficiency for stopped nucleons we use our data at the 
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Figure 9: The distribution of total energy deposition in calorimeters for candidates to 
e + e~ —> pp events with c.m. energy below threshold (left) and above threshold (right) 


energy point 950 MeV (run 2012), where part of antiprotons stop in the 
beam pipe and in the DC inner shell, and annihilate, N ann , while part of pp 
pairs pass the DC sensitive volume, and are identified as collinear events, 
N coi i, allowing to calculate the cross section. A fraction of stopped and 
annihilated antiprotons vs beam energy e s t opP ed was obtained from simula¬ 
tion and is shown in Fig. CCD by circles. Using the measured beam energy ( 
E c .m. = 950.8 MeV see above ), we obtain chopped = 0.5 ±0.1 with the corre¬ 
sponding e co u = 0.20 ±0.04 for the detected pp collinear events (Fig. [6]). (The 
uncertainties of e st0 pped was obtained from simulation with different thickness 
of vacuum pipe, the thickness is 0.50 ± 0.05 mm.) At this energy point we 
calculate a “visible” cross section for the collinear events, <r vis = V co ii/(e co nL), 
and assuming the same production cross section, we determine a detection 
efficiency for the annihilated antiprotons, e ann : 


AC 


^ann * ^stopped 


^ann ^coll-^ ^ann ^coll 


L * 0"vis 


L ■ N col i 


N, 


-± 


coll 




^ann ^coll 


(44.8 ± 9.2) • (0.20 ± 0.04) 
Acoiiestopped = (164.6 ±13.0) -(0.5 ±0.1) 


= 0.112 ±0.033 
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Figure 10: Numbers of e + e — > pp events with annihilation in the beam pipe and in the 
DC inner shell. The vertical line shows the proton-antiproton production threshold. 


The detection efficiency for annihilated antiprotons is found to be e ann = 
0.112 ±0.033, and because antiprotons annihilate at rest, it does not depend 
on the beam energy. In the 2011 run the real beam energy at 950 MeV 
is a little below the threshold of DC penetration, and we can not use the 
procedure. At other energy points the cross section can be calculated as: 


^Born — 


N - 

1 y pp 


i^ann ^stopped (1 C) 


The obtained values are listed in Table [21 


4. The \Ge/Gm\ ratio 

The e + e _ —* pp cross section depends on the proton/antiproton polar 
angle 6 as: 




iuy 2 l3C 

2s 


ml 


|G m (s)| 2 (1 ± cos 2 9 ) H - -\G E (s)\ 2 sin 2 6 


(4) 


and the \G E /G 
angle distribution. 


m | ratio can be extracted from the experimental polar 

AlKI ^ 

The ratio depends on energy, also —-^IG^s)! 2 is not 
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Figure 11: Data from simulation for the fraction of antiprotons annihilating at rest in the 
beam pipe and in the DC inner shell (dots), and those, penetrating deeper than 15 cm in 
DC (stars). 


equal 1, but in spite of this fact we combine the angular distributions from 
all c.m. energy points in the 1920-2000 MeV interval, because of insufficient 
statistics. We use the procedure described in Sec. 13. II to obtain the data-MC 
correction for the angular dependence of the detection efficiency. 

These corrections are shown in Fig. [T2] for the two experimental runs. We 
fit the corrected experimental pp polar angle distribution shown in Fig. [[3] 
by points with a sum of two functions 


sim , 
b G E =0 + 


I Gi 


\G 


M\ 


rpsmi 

|2"' r G , M=0’ 


(5) 


where Fq£ =q and Fq™ =0 are contributions to the angular distribution 
obtained from simulation with Ge = 0 and with Gm = 0, respectively. The 
numbers of simulated events for Gm = 0 and Ge = 0 are normalized to the 
integral of the (1 + cos 2 6) and sin 2 9 functions, respectively. The fit yields 
\G e /G m \ = 1.49 ±0.23. 

We estimate a systematic error on this value as 20%, mostly coming from 
the large statistical errors of the angular correction to the efficiency of FigJI?l 
and also strong dependence of the corrections on F c m .. 
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Figure 12: The efficiency correction 

(e®g|/e“™) 2 vs antiproton polar angle for the 
two experimental runs. 



Cose 


Figure 13: MC-simulated proton polar angle 
distributions with Ge= 0 (dashed), Gm= 0 
(dashed-dotted), and fit of experimental 
data (solid line). 
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Figure 14: The |Ge/Gm| ratio found in 
this work in comparison with other ex¬ 
periments. 


Figure 15: The proton effective form fac¬ 
tor measured in this work and in the 
BaBar experiment 0] . 


A comparison of the measured \Ge/Gm\ value with other experiments is 
shown in Fig. [TH 

5. Systematic errors 

5.1. Systematic errors; pp are detected 

The main sources of the systematic errors in this energy range are: 

( £ sim \ 2 

-Kfr ratio accuracy - 5%; 

£ reg J 

- variation of the \Ge/Gm\ ratio within error bars leads to 3% changes in 
the detection efficiency; 

- selection criteria - 2%; 
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- luminosity determination - 1% ,15]; 

- radiative corrections - 1% 13]; 

- beam energy determination accuracy - less than 0.5% above 955 MeV 
and 2% at energy point 950 MeV (run 2012); 

- uncertainty of the beam pipe thickness - 2% above 955 MeV and 10% 
at 950 MeV (run 2012); 

Combining above numbers, we estimate a total systematic error as 6% for 
beam energies above 955 MeV and 12% for the energy point 950 MeV(run 
2012 ). 


5.1.1. Systematic errors; p annihilates 

Main sources of systematic errors in this energy range are: 

- accuracy of antiproton detection efficiency - 32%, including uncertainty 
on the beam energy determination and uncertainty in the beam pipe thick¬ 
ness ; 

- \Ge/Gm\ accuracy gives 8% errors in the number of antiprotons stopped 
in the beam pipe and the DC inner shell. This contribution is considered in 
antiproton detection efficiency; 

- luminosity determination - 1 %; 

- radiative correction - 1 %; 

- selection criteria - 2 %; 

The total systematic errors are 33%. 

The measured cross sections are shown in Fig.[15]and Fig. [T6l respectively. 
The cross sections are listed in Table |T] and Table [2] 


6. CONCLUSION 

The e + e _ —> pp cross section has been measured using a data sample 
of 6.8 pb _1 collected in the center-of-mass energy range from pp threshold to 2 
GeV. Results agree with the previous BaBar experiment and have comparable 
or better statistical and systematic errors. The value of the ratio \Ge/Gm\ = 
1.49 ± 0.23 ± 0.30 has been found in the energy range from 1.92 to 2 GeV, 
and is in agreement with BaBar data. The expected tenfold increase in the 
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Table 1: The c.m. energy, beam energy shift, luminosity, number of selected e + e _ —> pp 
events, detection efficiency, radiative correction, and cross section with statistical and 
systematic errors. The data for collinear type events. 


£ c . m .,MeV 

«n,MeV 

L, nb 1 

N pP 

e 

(1-5) 

a, nb 

1900 (2012) 

0.8±0.1 

900.0 

164±13 

0.2 ±0.04 

0.75 

1.2±0.26±0.14 

1920 

3.3T0.2 

566.9 

251±16 

0.631 

0.81 

0.87±0.05±0.05 

1925 

0.5T0.3 

590.8 

280±17 

0.638 

0.82 

0.90±0.05±0.05 

1940 

2.4T0.4 

993.8 

488±22 

0.669 

0.85 

0.87±0.05±0.05 

1950 

1.2T0.3 

451.0 

238±16 

0.692 

0.86 

0.89±0.06±0.05 

1960 

3.0T0.3 

692.2 

397±20 

0.685 

0.87 

0.96±0.06±0.06 

1975 

1.3T0.3 

506.6 

283±17 

0.708 

0.88 

0.90±0.05±0.05 

1980 

3.6T0.5 

600.6 

356±19 

0.693 

0.88 

0.98±0.05±0.06 

2000 

2.3T0.4 

478.0 

284±17 

0.708 

0.88 

0.95±0.06±0.06 

Total 


5770. 

2741±52 

— 

— 

— 


Table 2: The c.m. energy, luminosity, number of signal events, fraction of antiprotons 
stopped in beam pipe and DC inner shell, efficiency, cross section with statistical and 
systematic errors, for annihilation events. 


^c.m.jMeV 

z^shitt 

■^beam 


N pP 

(1-5) 

^ stopped Cann 

<7, nb 

1890 

0.0 ±0.5 

527.1 

79.4±11 

0.69 

0.110 

1.98 ±0.27 ±0.66 

1900 (2011) 

0.0 ±0.5 

498.5 

41.3±8.5 

0.75 

0.067 

1.65 ±0.34 ±0.54 

1900 (2012) 

0.8 ±0.1 

900.0 

44.8±9.2 

0.75 

0.055 

1.21 ±0.25 ±0.40 
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Figure 16: The e + e —> pp cross section measured in this work (circles) in comparison 
with the BaBar [7) data (triangles). Only statistical errors are shown. 
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luminosity of VEPP-2000 will allow to measure the proton form factor and 
\Ge/Gm\ ratio with much better accuracy. 
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